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The nonequilibrium hydrodynamic correlations of a Multiparticle-Collision-Dynamics (MFC) fluid 
in shear flow are studied by analytical calculations and simulations. The Navier-Stokes equations 
for a MFC fluid are linearized about the shear flow and the hydrodynamic modes are evaluated as 
an expansion in the momentum vector. The shear-rate dependence and anisotropy of the transverse 
and longitudinal velocity correlations are analyzed. We demonstrate that hydrodynamic correlations 
in shear flow are anisotropic, specifically, the two transverse modes are no longer identical. In 
addition, our simulations reveal the directional dependence of the frequency and attenuation of the 
longitudinal velocity correlation function. Furthermore, the velocity autocorrelation functions of a 
tagged fluid particle in shear flow are determined. The simulations results for various hydrodynamic 
correlations agree very well with the theoretical predictions. 

PACS numbers: 


I. INTRODUCTION 

The thermodynamics of systems far from equilib¬ 
rium has drawn growing interest in the last couple of 
decades [I|. Several non-equilibrium relations, collec¬ 
tively called fluctuation relations, have been derived for 
transient and steady non-equilibrium states. These re¬ 
lations have been verified using exactly solvable mod¬ 
els and numerical simulations (see Ref. [l| and refer¬ 
ences therein). An interesting class of non-equilibrium 
systems is fluids under external fields such as shear 
flow and/or a temperature gradient 0- Gonsiderable 
progress has been achieved in understandiM these sys¬ 
tems using hydrodynamics calculations [3l-l8|, numerical 
simulations (9l-[iT|. and experiments [l^. For instance, 
the fluctuation relation for entropy production has been 
verified in numerical simulations of simple fluids under 
shear flow [l^, [l^ . Apart from satisfying fluctuation re¬ 
lations, non-equilibrium fluids show several interesting 
features which are absent in equilibrium. In particu¬ 
lar, non-equilibrium hydrodynamic correlations in steady 
states are long-ranged even for fluids far from critical 
points d, 0, S 0, Il3| . In addition, these correlations 
are anisotropic, in contrast to equilibrium correlations in 
simple fluids. A consequence of the long-range nature of 
the correlations is the non-intensivity of pressure fluctu¬ 
ations d]. 

Computer simulations are extremely valuable to 
study nonequilibrium phenomena. In particular, re¬ 
cently developed mesoscale hydrodynamic simulations 
approaches, such as lattice Boltzmann [iM3> dissipa¬ 
tive particle dynamics (DPP) or multiparticle 

collision dynamics (MPC) permit to cover large 

length and long time scales, and a wide range of external 
parameters such as shear rates and temperature gradi¬ 
ents. All the approaches are essentially alternative ways 
of solving the Navier-Stokes equations for the fluid dy¬ 
namics. Common to them is a simplified, coarse-grained 


description of the fluid degrees of freedom while main¬ 
taining the essential microscopic physics on the length 
scales of interest [l^. By now, the MPC method has suc¬ 
cessfully been applied in a broad range of equilibrium and 
nonequilibrium simulations of soft matter systems (see, 
e.g.. Ref. [ill and references therein). In particular, the 
hydrodynamic correlations of the MPC fluid have been 
determined and it has been shown that they agree with 
the solutions of the fluctuating Landau-Lifshitz Navier- 
Stokes equations [il|. Moreover, the hydrodynamic cor¬ 
relations of embedded colloids and polymers 

have been calculated. Even more, MPC simulations have 
been successfully applied to verify the fluctuation rela¬ 
tion for entropy production in shear flows EH- So far 
however, an analysis of nonequilibrium correlation func¬ 
tions of a MPC fluid and a comparison with theoretical 
approaches is missing. 

In this paper, we fill this gap and determine analyti¬ 
cally and by MPC simulations the time-correlation func¬ 
tions of hydrodynamic variables of a simple isothermal 
fluid under shear flow. We first derive analytical ex¬ 
pressions for the respective correlations by linearizing 
the Navier-Stokes equations. To this end, we follow the 
methods employed in Refs. [1, Q, where adiabatic or 
granular fluid are considered. Here, the isothermal ap¬ 
proach is simpler, because energy is no longer a conserved 
quantity. We restrict ourselves to moderate shear rates 
for which the coupling between hydrodynamic modes can 
be ignored SO- Exploiting the MPC method, we then 
perform shear flow simulations and calculate the respec¬ 
tive hydrodynamic correlation functions. The primary 
effect of shear is the anisotropy of the hydrodynamic cor¬ 
relation functions, as already predicted in Refs. iO- 
The frequency and attenuation of the longitudinal modes 
become directional and shear rate dependent. In addi¬ 
tion, the degeneracy of the two transverse modes, present 
at equilibrium, is removed. The anisotropy of the longi¬ 
tudinal and transverse velocity autocorrelations is also 
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manifested in the anisotropy of the velocity autocorrela¬ 
tions of tagged MFCs particles. Moreover, the correla¬ 
tion functions show a faster decay than the equilibrium 
correlations at long times. By comparison, we find excel¬ 
lent agreement between the theoretical predictions and 
the MFC simulation results. 

The article is organized as follows. The theoretical ex¬ 
pressions for the velocity correlation functions are derived 
in Sec. im Section m presents simulation results and a 
comparison with the theoretical predictions. Our results 
and findings are summarized in Sec. IIVI More details of 
the calculations are presented in the Appendices. 


II. THEORY 


A. Linearised Navier-Stokes equations under shear 


The Navier-Stokes equations of an isothermal MFC 
fluid are given by 


dp 


= -V • (pu) 


( 1 ) 




u V 


U =-Vp-f -f (V • u) . (2) 


They account for mass and momentum conservation, 
where p(x, t) is the mass density, u(x, t) the fluid ve¬ 
locity field, and p(x, t) the pressure field at the posi¬ 
tion X at time t. The shear viscosity is denoted as 
77 . The Navier-Stokes equations are adopted to a non- 
angular-momentum-conserving MFC fluid, hence, the ki¬ 
netic contribution 77 ^ of the shear viscosity appears in the 
last term in the rhs of Eq. rather than the viscosity 
Tj itself [ 2 ^. In addition, we omit the fluctuating part 
of the stress tensor ii in Eq. ©■ The equations are 
then linearised by setting p = po + F = Fo + Sp, and 
u = Uo-l-(5u, where uoa = lapxp, with the shear-rate ten¬ 
sor 7 o,p and a, /3 G {x, y, z}. We choose the x- and p-axis 
of the Cartesian coordinate system as the flow and the 
gradient direction, respectively, such that 'jap = jSaxSpy, 
where 7 is the shear rate. We use the summation conven¬ 
tion for Greek indices unless otherwise stated. Equations 
0 and ([2]) can then be written as 


d 


d 


Po 




d d 


6p = — poV • (5u 


(3) 


d 


SUa = - Po^apSup - 


-|-77V^(5Mo 


77 ^ d 
3 dxo 


(V • 5u) . (4) 


Here, we have neglected second order terms in the fluc¬ 
tuations. We eliminate 5p with the ideal gas equation of 
state, 5p = c^6p, where ct is the isothermal velocity of 
sound. By rescaling the velocity and density according 
to du = i5u/ct and 5p = (5p/po, Eqs. ([2|) and (|3]) can be 


written in momentum space as 

5p =zcTk • (5u 


d , d 


(5) 


d , d 


dt 


Sua = - + icrkaSp 


— VK SUoL - —kakjS Suj3 

o 


( 6 ) 


with the kinematic viscosities 1 / = p/po, = 77 ^/po. The 
variables with a tilde are Eourier-transformed variables 
according to the definition 


f(k) 


d3xe**"'’"f(x). 


(7) 


We now write the above equations in terms of the lon¬ 
gitudinal and transverse component of the velocity field. 
Let (5u = where 

and are three orthogonal unit vectors. Here, is 
chosen along the propagation direction of k, so that 
is the longitudinal, and and are the transverse 
component of the velocity field. By introducing the vec¬ 
tor z = {Sp, , the Navier-Stokes equa¬ 

tions can be written as 




ikx 


dk, 


Vi 


i + a = Q. 


( 8 ) 


The explicit form of the matrix C for the choice Q 


= k/|k| 

= [y - eWe(i)]/fc_L (9) 

e(3) = X e(2) 


of the unit vectors is given in Appendix O Here, y is the 
unit vector along the p-axis in the Cartesian coordinate 
system and k±_ = -f where k = |k|. The so¬ 

lution to the above equation can be written as the linear 
combination 


z (k, 0 = X] 

of the eigenvectors which satisfy the eigenvalue 

equation 

[-7fc.^+£]|(*Hk) = A^^*nk). (11) 

Oky 

Let J 7 (®)(k) be the corresponding left eigenvectors such 
that 






( 12 ) 
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The left and right eigenvectors and the eigenvalues can be 
calculated using perturbation theory Q and are given in 
Appendix. IbI Inserting z(k, t) from Eq. (|TTH) into Eq. (|8]) 
and using Eq. (ED together with the orthogonality con¬ 
dition in Eq. dE, we obtain 

The solution of the above equation is given by 

a(*^(k,t) = a(*^(k(-t),0)exp drAi(k(-T))^ , 

(14) 

where the time dependent k vector is defined as k(t) = 
{kx,ky — jtkx,kz). Using Eqs. (fTUl) and (ITTl) and the 

relation a(*)(k, 0) = J2t=i 0); 8®^ 

4 

ZiiKt) ='^Gi3{Kt)zj{H-t),0), (15) 

1=1 


nonzero for 7 ^ 0. We consider only small shear rates 
7 ^ vk^, for which the cross-correlations can be ne¬ 
glected. Hence, the correlation functions can be writ¬ 
ten as Cij{k,k',t) ~ (27r)^(5yd(k-t-k')Cii(k,<), where 
Cii{k,t) = (zi(k,0)zi(-k,0)) Gii{k{t),t). Using the ex¬ 
plicit expressions for the propagators Gii{k{t), t), the cor¬ 
relation functions can be written as 


Cii(k,t) 


Crp k J 


C22(k,<) 


C33(k,t) 


-|■C'll(k, t), 

Po 

ksT f \ 

Po \k{t)) 




C44(k,t) = 

Po 

where 4>{k,t) and x(k, t) are given by 


cos [ct<(> (k, t)] , 
( 20 ) 

( 21 ) 

( 22 ) 

(23) 


where the propagator Gy (k, t) is defined as 


Gy (k, t)=J2 (k)i?® (k(-t)) exp 

Z=1 


drAi(k(-T)) 


ae) 

In order to compare with the simulations, it is convenient 
to rewrite Eq. (TTSl) by setting k = k(t). We then get 


cl){k,t) = 


2'jkx L 


— In 


[kyk — ky{t)k{t)] 
ky{t) + k{t) 


ky + k 


X(k, t) = kH - -^kxkyt^ + ^i'^klt^, 


(24) 

(25) 


Zi{k{t),t) =^Gy (k(t),t) %(k,0), (17) 

1=1 


with 


Gy (k(t),t) = ^ef^(k(t))77]'^(k)exp (^- drAz(k(T))^ , 

(18) 


using 


fo dTA,(k(-T)) 


. k=k(t) 


= fg (irAi(k(r)). The ex¬ 
plicit form of Gij (k(t), t) can be obtained from the 
eigenvectors {^,17} and the eigenvalues A’s given in Ap- 
pendix[BJ Note that the solution given by Eq. (fT71) repre¬ 
sents the evolution of the hydrodynamic variables in the 
time-dependent reference frame in the k-space. 


B. Hydrodynamic correlation functions 

The correlations of the hydrodynamic variables are de¬ 
fined as Gy(k, k',t) = (zi(k(t), t)5j (k', 0)), and become 

with Eq. ED 

4 

Gy(k,k',t) = ^G,4k(t),t) (5,(k,0)5,(k',0)). (19) 

1=1 

The correlations (zi(k,0)zj(—k,0)) vanishes at equilib¬ 
rium, i.e., 7 = 0, for i ^ j. However, they are 


and k{t) = |k(t)|. Here, v = v + /?>, and the equi¬ 
librium relations (5i (k, 0) (—k, 0)) = poksTc^^ and 

(5i(k,0)2i(~k,0)) = pg^ksT for i = 2,3,4 have been 
employed. These expressions can be derived using fluc¬ 
tuating hydrodynamics for a MFC fluid [s^, however we 
do not present the derivations here. 

A few remarks on the correlation functions given by 
Eqs. (EOD-dMD are in order. In the limit 7 —> 0, we get 
(j){k, t) —7> kt and x(k, t) —>■ k'^t^ and therefore the correla¬ 
tion functions are reduced to the corresponding equilib¬ 
rium relations [2^, to 0{k^). In the absence of shear, 
the correlation functions for an isothermal MFC fluid can 
be obtained for all orders in k; the exact expressions for 
the velocity autocorrelations are provided in Ref. [^ . 
We also note that the expression for G33(k, t) remains 
exact for all shear rates within the order we are working 
at, even if the neglected equal-time correlations of the 
form (5i(k, 0)5j(—k, 0)) for i j are taken into account. 
By the same token, G44 (k, t) is exact for all shear rates 
for kz = 0. 


C. Velocity correlations in real space 


From Eq. ED, the velocity correlation function follows 
as 

((5u(k, t) ■ (5u(k', 0)) = (27r)^(5(k(—t) -|- k')G“(k, t), (26) 
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with the abbreviation 

4 

C“(k,t) = '£CM-t),t) e«(k) • e«(k(-t)) (27) 


i=2 


and by using (^^(k, t) ~ 0 for * j. The velocity auto¬ 
correlation is real space is then given by 


((5u(x, t) ■ du(0, 0)) = 


1 


(2^)^ 


d^k C“(k,t)e 


— ikx 


(28) 


The velocity autocorrelation function C{t) = (v(t) • v(0)) 
of a tagged particle of velocity v(t) can be obtained by 
setting v(f) = u(r,t), where r is the position of the 
tagged particle, and averaging over all its positions r. 
Hence, we obtain 


C{t) = 


1 


(2^) 


d^k C'“(k,t)(e-*‘^ '’) 


(29) 


with the definition (e®*^ '’) = / drP(r, and P(r, t) 

the distribution function of the position of the tagged 
particle. Using the Fourier representation of P(r, f), we 
get (e“*'^''’) = P(k,t). In shear flow, P(k, t) follows from 
the advective diffusion equation [s^ 

P{k,t) =-Dk^P{k,t), (30) 


d d 

^rC; 


dt 


dky 


where D is the diffusion coefficient. The solution of the 
equation is 

ct 


P(k, t) = P(k(—t), 0 ) exp U J dTk‘^{—T) 


Then, Eq. (|2^ yields 
1 


C{t) = 


( 27 r) 


Pk C'“(k,f)P(k(-t),0) 


(31) 


(32) 


X exp J dTk^{—T)^ . (33) 

By changing the integration variable from k to k(t) , and 
using the fact that the Jacobian of the transformation is 
unity, we get 

C{t) = - 3 f d^k C'^{k{t),t) exp (—D f dTk^{T) 

( 27 r) J \ Jo 

(34) 

by using P(k,0) = 1 (H. 

So far, we considered infinitely large systems. In com¬ 
puter simulations, however, finite-size systems are used 
with typically periodic boundary conditions. This leads 
to a discrete set k^ of wavevectors, with k^^n = STTric/P, 
where L is the length of the cubic simulation box of vol¬ 
ume V = L^, Ua S Z, and kn ^ 0. Hence, the correlation 
function becomes 

C^{^n{t),t)exp(-D J dr k^ (t)V 

(35) 



t/T 


FIG. 1: (Color online) Numerically evaluated transverse ve¬ 
locity correlation function along e^®^(k) (see Eqs. ([51 and 
(1231)). The lines (solid) correspond to shear rates yr = 
0.0, 0.001, 0.01, 0.1, 1.0 (right to left). In the main figure, the 
wavevector components are fcx = 27r/60 and ky = kz = 0 , 
and in the inset kx = ky = 27r/60 and kz = 0. 


The velocity autocorrelation function in shear flow is 
anisotropic. Therefore, we write the above equation in 
terms of the components corresponds to the three or¬ 
thogonal directions as 

-l 00 3 

E E^p(k„,t)e«(k„)e«(k„(t)), 

kn = —00 1 — 1 

(36) 

where j = Z -|- 1. Note that index a is not summed 
over. Since MFC is a particle-based mesoscale simula¬ 
tion method, the validity of the Navier-Stokes equation 
breaks down at the level of collision cells [ 2 ^. Therefore 
the k-values in the summation in Eq. (I36p are limited 
by a cut-off corresponding to the smallest hydrodynamic 
length scale. Alternative but similar approaches to evalu¬ 
ate the velocity autocorrelations of a tagged fluid particle 
can be found in Refs. S®- 


III. SIMULATIONS 

A. Multiparticle collision dynamics 

In the MF C ap proach, the fluid is represented by point- 
particles (^ . . Their time evolution proceeds in two 

independent steps, namely the streaming and collision. 
In the streaming step, the particles move ballistically, 
i.e., the particle positions are updated as 

Xi{t + h) = Xi{t) + hvi{t) , (37) 

where h is the collision-time step. Here, x^ denotes the 
position of particle i, its velocity, and i € {!,..., N}, 
with the total number of particles JV. In the collision 
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FIG. 2: (Color online) Theoretical and simulation results 
for the transverse velocity correlation functions along (k) 
(blue, dotted) and e^^^(k) (red, dashed) for shear rates yr = 
0.0, 0.005, 0.01 (top to bottom at t/r = 100). The two trans¬ 
verse components are identical for -yr = 0.0, and therefore 
only one of them (red, dashed) is presented. The solid lines 
(black) represent the theoretical results. In the main figure, 
the wavevector components are =2-k jL and ky = kz = 0, 
and in the inset k^ = ky = 21^ jL and kz = 0. 


step, the particles are grouped into cubic cells of length 
a, and a rotation of their relative velocities—with respect 
to the center-of-mass velocity of the particular cell—is 
performed. Hence, the new velocities are 


Vi(t + h)= Vcm{t) + Tl{a) [Vi(t) - Vcm(t)] (38) 


Here, Vcm{t) is the center-of-mass velocity of the cell that 
contains the particle i and TZ{a) is the rotation matrix, 
with the axis of rotation taken as a random unit vector. 
A random shift of the collision cell lattice is performed at 
every collision step to ensure Galilean invariance [2^, [s^ . 

We perform isothermal simulations, where tempera¬ 
ture is maintained by the cell-level Maxwell-Boltzmann- 
Scaling (MBS) approach, which has been shown to yield a 
canonical ensemble . The hydrodynamic fluctua¬ 

tions of the MFC fluid supplemented by the MBS method 
is known to be consistent with the linearized Navier- 
Stokes equation in equilibrium [23, HI]. Shear flow is 
implemented by Lees-Edwards boundary conditions [s^ . 
The time step of our simulations is chosen as h/r = 0.1, 
with the unit of time r = y/ma^/kBT, to ensure a large 
Schmidt number [s^, and the average number of parti¬ 
cles in a collision cell is set to 10. The numerical values 
of the transport coefficients for this choice of the sim¬ 
ulation parameters are v = O.STOa^/r, zz = O.SSTa^/r, 
D = O.OSlaVr, ct = l.Oa/r (Hj. 


B. Hydrodynamic correlations 

The density and velocity fields in k-space are defined 
as 

N 

p(k,t)=^e*‘^W-% (39) 

i=l 

N 

5u(k,t) = [vi - Uo(x,;)] (40) 

where uq(x) = yyx is the mean velocity field. 
Note that we use the time dependent k-vector 
k(t) = {kx,ky — jtkxTkz), so that the definitions of the 
hydrodynamic fields in momentum space are consistent 
with the Lees-Edwards boundary conditions. The trans¬ 
verse and longitudinal components of the velocity held 
are then dehned as 5u^®^(k(t), t) = e*^*^(k(t)) • (5u(k, t). 

A few notes on the calculation of autocorrelation func¬ 
tions in shear how implemented via Lees-Edwards bound¬ 
ary condition are in order. In equilibrium simulations, 
the origin of time is arbitrary, and therefore the moving- 
time-origin scheme [40l| for calculating time correlation 
functions can be employed to improve statistics and to 
avoid storing position and velocity coordinates of the par¬ 
ticles. However, in our simulations, the k-vector is taken 
as a function of time and the time origin is taken as 
the time at which the image of a particle in the inhnite 
periodic system is given by x' = x^ -|- L, where x^ is the 
position of the particle in the primary simulation box and 
L = L{nx, Uy, nzf" ■ Therefore, averages have to be taken 
only over the allowed time origins. In addition, in order 
to be consistent with the dehnition of the time-dependent 
k-vector, the position coordinate in the gradient direction 
has to be taken in the range \—Ly/2, Ly/2]. However, the 
usual moving-time-origin scheme can be employed in the 
evaluation of real-space time-correlation functions. 


1. Correlation functions in momentum space 

Figure [T] shows the numerically evaluated transverse 
velocity correlation function given by Eq. (USD. As is ev¬ 
ident from the theoretical expression, there are primarily 
two time-regimes for the decay of the correlations. For 
t <C 1 / 7 , the decay is dominated by the term linear in t in 
the exponential (see Eq. (1^51) ) and therefore is identical 
to the decay of the correlation function in equilibrium. 
However, for t » I/ 7 , the decay is dominated by the 
term proportional to and is characteristic of the shear 
flow. This term originates from the advection term in the 
Navier-Stokes equation and results in several interesting 
features such as faster decay with a power-law of 

the long-time tail in the velocity autocorrelation function 
of a tagged fluid particle and renormalization of the 
viscosity [^ . 
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FIG. 3: (Color online) Longitudinal velocity correlations for 
the shear rates 7 r = 0.0 (red, dashed) 0.005 (blue, dotted), 
and 0.01 (green, dotted-dashed). The solid lines (black) rep¬ 
resent theoretical results. The wavevector components are 
kx =2 -k/L and ky = kz = 0. 


In Fig. [5] we compare the transverse velocity corre¬ 
lations obtained from the simulations and the theoret¬ 
ical expressions. In contrast to equilibrium correlations, 
the autocorrelations of the two transverse components 
in shear flow are not identical. The transverse velocity 
component perpendicular (e^^^ direction) to the gradi¬ 
ent direction decays slower than the second component 
perpendicular direction) to the longitudinal direc¬ 
tion for long times (t ^ l/b)- Even though the distinc¬ 
tion between the two transverse components is apparent 
from our simulations, it may be ignored in deriving the 
long-time tail exponents for the velocity autocorrelation 
function of a tagged fluid particle The transverse 



FIG. 4: (Golor online) Longitudinal velocity correlations for 
7 T = 0.005, kx = ky = 27r/L (dotted-dashed), 7 r = 0.005, 
kx = —ky = 27r/L (dotted) and 'yr = 0.0, kx = ky = 27r/I/ 
(dashed), kz = 0 for all the curves. The solid (back) lines 
represent the theoretical results. 


correlations, as mentioned in the previous paragraph, de¬ 
cay similar to that in equilibrium for t <C I /7 and faster 
for t :$> 1 / 7 . For t Ri I/ 7 , the decay depends on the 
direction of the k vector, i.e., on the relative sign of kx 
and ky. For sgn(fca;) = sgn(fcy), both the transverse cor¬ 
relations decay slower than the equilibrium correlations, 
and faster otherwise (see insets of Figs. [1] and For 
ky = 0, the transverse correlation functions decay faster 
than the equilibrium correlations at all times. 

The longitudinal velocity correlation function corre¬ 
sponds to the sound propagation in the fluid. There are 
two effects of shear flow on the propagation of sound in an 
isothermal fluid - the modification of the sound damping 
factor and the change in the sound frequency/velocity 
(Doppler effect), both of which depend on the shear rate 
and the direction of propagation. Figure [3] shows the 
variation of longitudinal velocity correlations in the flow 
direction {ky = kz = 0) for different shear rates. The 
change in the frequency and the faster attenuation with 
increasing shear rate is well demonstrated. Figure 0] dis¬ 
plays the anisotropy of the sound propagation. The fre¬ 
quency decreases when the sound propagation is in the 
direction along the flow, and increases in the direction 
against the flow. The direction dependence of the atten¬ 
uation of longitudinal velocity correlations is the same 
as that of the transverse velocity correlations. The au¬ 
tocorrelation function of the density fluctuations shows 
an identical behavior as the longitudinal velocity corre¬ 
lations, and therefore we do present the results here. 


2. Long-time behavior of velocity correlations 

Figure [5] shows velocity autocorrelation function of a 
tagged particle. Note that we consider the thermal ve¬ 
locity of the particle, i.e., the velocity with respect to 
the mean flow velocity. Evidently, the correlations in the 
three orthogonal directions are not identical. We find ex¬ 
cellent agreement between theory and simulation results 
for long times. The deviations at short times are caused 
on the one hand by the fact that the theoretical hydro- 
dynamic correlations are only accurate to 0{k‘^). On the 
other hand, partition of the MFC fluid in collision cells 
leads to a break-down of hydrodynamics at short times 
and length scales below the collision-cell size [ 2 ^. How¬ 
ever, the long-time behavior is determined by small k 
values, i.e., large length scales, which are correctly repro¬ 
duced in the simulations. 


IV. SUMMARY AND CONCLUSIONS 

We have studied the nonequilibrium hydrodynamic 
time correlations of an isothermal MFC fluid under shear 
flow. We find good agreement between simulation re¬ 
sults and theoretical predictions based on the linearized 
Navier-Stokes equations for moderate shear rates. We 
confirm that hydrodynamic correlations in shear flow are 






FIG. 5: (Color online) Velocity autocorrelation functions of a tagged particle along the various spatial directions. The shear 
rate is 7 r = 0.01. Simulation results are represented by open circles and the theoretical prediction by solid lines. 


anisotropic, in agreement with previous studies iii. 
Specifically and contrast to equilibrium correlations, the 
time correlations of the two transverse modes in Fourier 
space are no longer identical. In addition, our simulations 
reveal the directional dependence of the frequency and 
attenuation of the longitudinal velocity correlation func¬ 
tion. As a consequence, the velocity autocorrelation of 
a tracer fiuid particle (MFC particle) is also anisotropic. 
The agrement between analytical calculations and simu¬ 
lations confirms that MFC is a suitable approach to study 
hydrodynamic properties of simple fluids under nonequi¬ 
librium conditions. 


where £ = —ikCi + k'^C 2 + 7 /I 3 , with 

/ 0 CT 0 0\ 

^ _ ct 0 0 0 I 

1 “ 0 0 0 0’ 

Vo 0 0 0/ 


/O 0 0 0\ 

0 z> 0 0 I 

0 0 Z/ 0 ’ 

Vo 0 0 z// 


(A2) 


Our studies are restricted to moderate shear rates, 
where equal-time correlations of the hydrodynamic vari¬ 
ables can be approximated by the corresponding equilib¬ 
rium values. For high shear rates, we observe significant 
deviations of the simulation results from the theoreti¬ 
cal expressions. The deviations increase with the shear 
rate. In order to theoretically evaluate the equal-time 
and autocorrelation functions for high shear rates, the 
fluctuating part of the stress tensor has to be included 
in the Navier-Stokes equations Q, which we omitted. In 
addition, it is also necessary to take into account the den¬ 
sity dependence of the viscosity in linearising the Navier- 
Stokes equation. These issues will be addressed in future 
publications. 



0 

0 


0 

Til 

ri2 

£13 

0 

r2i 

£22 

£23 

Vo 

Tai 

£32 

£33/ 


where z> = z^ -|- z^^/3, and the matrix F is defined as 


r\ 

b/T ■ ■ - — pOl^ lie _p(t) 

ji. Iml’-'i P-n 

(A3) 

For the particular choice of the unit vectors eO) 
in Eqs. (jH]), the matrix F takes the form 

as given 

/ kxky/k'^ 2kxk±/k'^ 0 \ 


F = 1 -kx/k± -k^ky/P 0 j , 

(A4) 

\—kykz/kk± —fcz/fc 0 / 



where k'j_ = k^ + k^. 


Appendix A: The hydrodynamic matrix 


The evolution of the hydrodynamic variables are given 
by 






z -b £z = 0 , 


(Al) 


Appendix B: Eigenvalues and eigenvectors of £ 

The eigenvalue equation Eq. (HU can be solved pertur- 
batively by expanding and Am in powers of k 

^{m) ^ 

Xjn — k Xm,l ■■ (^ 1 ) 












The solution to the order 0{k'^) is given by 
Ai = -icrk + i [vk'^ + '^k^ky/k'^) , 

A2 = +icTk + i [vk'^ + '^kooky/k'^) , 

A 3 = vk^ — 'jkxky/k'^, A 4 = vk^, (^2) 

= ^(1,1,0,or, = ^(1,-1,0,0)^ 

|(3) = (0,0,l,Mf, r") = (0,0,0, If , 

where 

M(k.) = — ^ ^ arctan f. (B3) 

kxkj_ Y kj_ J 
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